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Abstract: If one wants to explore the properties of a dynamical system systematically one has to 
be able to track equilibria and periodic orbits regardless of their stability. If the dynamical system 
is a controllable experiment then one approach is a combination of classical feedback control and 
Newton iterations. Mechanical experiments on a paramctrically excited pendulum have recently 
shown the practical feasibility of a simplified version of this algorithm: a combination of time- 
delayed feedback control (as proposed by Pyragas) and a Newton iteration on a low-dimensional 
system of equations. We show that both parts of the algorithm are uniformly stable near the 
saddle-node bifurcation: the experiment with time-delayed feedback control has uniformly stable 
periodic orbits, and the two-dimensional nonlinear system which has to be solved to make the 
control non-invasive has a well-conditioned Jacobian. 
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1. INTRODUCTION 



One way to explore a nonlinear dynamical system in a 
systematical fashion is bifurcation analysis by continuation: 
one starts in a parameter region where one knows that 
a simple attractor exists (say, a stable periodic orbit) 
and then varies a system parameter p, checking at which 
parameter values the periodic orbit loses its stability or 
"disappears" . At these special parameter values the periodic 
orbit undergoes a bifurcation, and other invariant objects 
(equilibria, periodic orbits, tori) may branch off. Thus, by 
systematically tracking equilibria, periodic orbits and their 
bifurcations in the parameter space one can (to a good 
extent) classify the qualitative behaviour of the dynamical 
system. 

In this paper we analyse the stability of the algorithm for 
pseudo-arclength continuation in experiments introduced 
in [Sieber et al. (2008)] that considered periodic rotations 
of a vertically excited pendulum near a saddle-node 
bifurcation. For this example experiment we demonstrate 
(using simulations) that all parts of the algorithm converge 
uniformly. In section 3 we recap how one embeds a 
Newton iteration into pseudo-arclength continuation to 
study how periodic orbits depend on system parameters. We 
also briefly explain how this continuation is implemented 
for the excited pendulum using time-delayed feedback 
control (TDFC; Pyragas (1992)) for the periodic part 
of the problem and a Newton iteration for only two 
scalar variables. We show that, with this approach, the 
system with TDFC is uniformly stable near the saddle-node 
bifurcation (which is not the case for the classical TDFC), 



and that the Jacobian used in the Newton iteration is 
uniformly well-conditioned. 

2. BACKGROUND ON RELATED METHODS 

If a dynamical system is given in the form of a low- 
dimensional ordinary differential equation (or discrete map) 
one can apply specialized algorithms based on Newton 
iterations embedded into pseudo-arclength continuation 
(see Section 3 for an explanation and an example), which 
are available as software packages, for example, Auto 
or Matcont (see textbooks [Doedel (2007); Kuznetsov 
(2004)] for a detailed introduction). These algorithms 
have been successfully extended to problems where the 
models are delay differential equations (Dde-Biftool: 
Engelborghs et al. (2001); Pddecont: Szalai et al. (2006)), 
dissipative partial differential equations (LoCA: Salinger 
et al. (2002); Lust et al. (1998)), or high-dimensional 
systems (such as stochastic Monte Carlo simulations) with 
'essentially low-dimensional' dynamics (Kevrekidis et al. 
(2004)). One advantage of algorithms based on Newton 
iterations is that they work independent of the dynamical 
stability of the state they track. Hence, they are also able 
to track unstable periodic motions or equilibria. 

In contrast to the situation where one investigates a model, 
bifurcation analysis in experiments is typically done by 
parameter studies: one gradually varies a system parameter 
and observes the transients. When one observes a slow- 
down of the transients or a sudden 'jump' of the output it 
is likely that one has encountered a bifurcation (and its type 
can sometimes be deduced from the transient behavior). 
Using an electronic implementation of a Duffing oscillator, 



Langer and Parlitz (2002) show how one can automate this 
approach to trace out bifurcations. Anderson et al. (1999) 
for an eloctro chemical system and De Feo and Maggio 
(2003) for a Colpitts oscillator have demonstrated the use 
of Newton iterations and pseudo-arclength continuation in 
experiments (also detecting or continuing bifurcations). 
Both studies ran a system identification procedure in 
parallel to the experiment, determined the steady states 
(fixed points or periodic orbits) of the identified model, 
and used feedback control to drive the experiment toward 
the identified steady state. The accuracy of the results 
using this approach is limited not only by the measurement 
accuracy and the tolerances set in the Newton iteration 
but also by the accuracy of the system identification. This 
is a severe handicap because system identification is an 
inverse and, thus, ill-posed problem. Anderson et al. (1999) 
demonstrated their approach for simulations of the chemical 
system only, actual experiments are still outstanding. 

An alternative approach to continuation of steady states in 
experiments is via the use of feedback control mechanisms 
that are automatically non-invasive: washout filters [Abed 
et al. (1994)] and time-delayed feedback [TDFC; Pyragas 
(1992); Socolar et al. (1994); Kim et al. (2001)]. Whereas 
classical feedback control compares the output of a 
dynamical system to a given reference signal <f> and feeds a 
(typically linear) combination of the difference 4> — 4> back 
into the experiment, washout-filtered feedback control and 
TDFC do not require a given reference signal. Washout 
filtered feedback picks 4> as the solution of 

j> = R-(<l>-4>) 

where R is a stable matrix. TDFC picks </> using the 
recursion 

4>(t) = (1 - R) ■ 4>{t - T) + R ■ <P{t - T) (1) 
where < R < 1. Thus, whenever a dynamical system 
subject to washout-filtered feedback control has a stable 
fixed point this fixed point is also a (possibly unstable) 
fixed point of the uncontrolled dynamical system (the same 
is true for periodic orbits of period T using TDFC). This 
means that bifurcation analysis could in principle be based 
on these non- invasive feedback control techniques. One 
difficulty encountered, and intensively discussed in the 
case of TDFC [Nakajima and Ueda (1998); Just et al. 
(1997); Hovel and Scholl (2005); Fiedler et al. (2007)], is 
finding reasonable conditions which guarantee that the non- 
invasive feedback is actually able to stabilize a steady state 
of a dynamical system. Typically, even if one has designed a 
feedback control loop that is able to stabilize a steady state 
0* when one inserts </>» as the feedback reference signal 
(that is, </> = </>*) there is no guarantee that 4>* is also stable 
if we replace the reference signal using the recursion (1) 
(which would correspond to TDFC). 

3. PSEUDO-ARCLENGTH CONTINUATION - 
VERTICALLY EXCITED PENDULUM EXAMPLE 

If one shakes the pivot of a pendulum up and down 
harmonically with a frequency to and an amplitude p then 
the pendulum can show stable rotations for amplitudes 
larger than a certain critical amplitude po- A simple model 
for the mechanical pendulum is 

ml 2 9 + W + ml[g + w 2 psin(wt)] sin(6>) = 0, (2) 
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Fig. 1. Bifurcation diagram for rotations of a paramet- 
rically excited pendulum at low forcing amplitude 
(experimental results from [Sieber et al. (2008)]). The 
y-axis shows the average phase (j>o as found by (14). 
Square marker: start of continuation; full circles: sta- 
ble part of the branch of rotations; large circle: saddle- 
node bifurcation occurs in uncontrolled pendulum; 
hollow circles: unstable part of family of rotations; 
diamond: loss of control. 

where 9 is the angle, m is the effective mass of the pendulum, 
I is its effective length, b is a viscous damping coefficient 
(typically small), g is the acceleration due to gravity, to is 
the excitation frequency in rad/s, and p is the exciation 
amplitude in m. For the parameters used in Fig. 1 the 
stable rotation at excitation amplitude p « 2 cm (square 
marker) is relatively easy to find by trial and error in a 
simulation or an experiment because it has a large basin 
of attraction. 

If one wants to find the minimal amplitude po that supports 
rotation one would in a simulation or experiment decrease 
the excitation amplitude in small steps, always waiting 
until transients decay after each parameter change. In this 
way one finds that the stable periodic orbit "disappears" 
at (or, rather, slightly above) the minimal amplitude po 
corresponding to the saddle-node in Fig. 1 as transients 
escape to another stable attractor, for example, the hanging- 
down state. In a study of a model, such as equation (2), 
the alternative to this vary-and-wait approach is a direct 
solution of (2) in rotating coordinates <f> = 6 — cut: 

mPfi + btft + bLd + mllg + ^psin^tysm^ + ujt) = 0, (3) 
with periodic boundary conditions 

0(T)-#O)=O, 0(0) - 0(T) = 0, (4) 
where T = 2it /lj is the (known) period of the rotation. The 
periodic boundary value problem (3), (4) is nonlinear and 
is typically solved with a Newton iteration. The advantage 
of this direct approach is that one can find periodic 
rotations independent of their dynamical stability: the 
periodic rotation undergoes a saddle-node bifurcation at the 
parameter value po such that the scenario looks as shown in 
Fig. 1. The Newton iterations for the nonlinear boundary 



value problem (3), (4) finds both, dynamically stable and 
unstable rotations in Fig. 1. Two difficulties for Newton 
iterations are: (i) it converges only locally, that is, a good 
initial guess is necessary; and (ii) the nonlinear problem 
is singular at the saddle-node in Fig. 1. Both problems 
can be overcome by embedding the Newton iteration into 
a pseudo-arclength continuation (see textbooks [Doedel 
(2007); Kuznetsov (2004)]): one treats the bifurcation 
parameter p also as a variable, such that the solutions 
(p, </>(•)) of (3), (4) form a curve C in the space of all possible 
functions and parameters. Furthermore, one extends the 
nonlinear boundary value problem by the (scalar) pseudo- 
arclength condition 



&an (*)[<£(*) " 0old] + <W*M*) - 0old(*)] dt- 



Tj 

+ Ptan[p-Pold] = h. (5) 

In (5) (poid) < / , oid(')) i s the previously found point on the 
curve C, (ptan, <Kan( - )) ^ s the unit tangent vector to the 
curve C in this previous point and h (a small quantity) is the 
approximate distance between (p id, &>id( - )) and the desired 
solution (p, </>(■)) of (3)-(5). Figure 1 shows a projection 
of this curve C. Using the pseudo-arclength extension (5) 
the nonlinear boundary value problem is uniformly well- 
conditioned along the whole curve C including the vicinity 
of the saddle- node at p — po - Pseudo-arclength continuation 
is a useful (and, by now, well established) tool in the 
numerical analysis of bifurcations because it allows one to 
follow unstable parts of branches of periodic orbits and also 
direct continuation of bifurcations (such as the saddle-node) 
in more than one parameter. In this way one can construct 
maps in the parameter space that help to classify for any 
given system its possible equilibria and periodic orbits, 
and their bifurcations [Krauskopf et al. (2007); Kuznetsov 
(2004)]. 

4. CONTINUATION IN EXPERIMENTS 

Extension of the pseudo-arclength continuation to experi- 
ments would give experimenters the opportunity to study 
bifurcations in much greater detail. For example, for the 
vertically excited pendulum with rotations as shown in 
Fig. 1 a parameter study that simply observes transients 
loses the stable periodic orbit already at amplitudes p signif- 
icantly above po (i n the preparatory studies for [Sieber et al. 
(2008)] at p = 0.8 cm for id = 3 Hz) due to disturbances 
or insufficiently small parameter steps. Thus, it would be 
difficult to establish that the loss of the stable periodic orbit 
is indeed due to a saddle-node bifurcation. Moreover, the 
average phase of the rotation with respect to the excitation 
changes dramatically within a tiny parameter range: a 
stably rotating pendulum points (nearly) upward whenever 
the pivot excitation reaches its maximum (at t — tt/(2lo)) 
whereas close to the saddle-node bifurcation the pendulum 
is nearly horizontal at time t = 2tt/uj. (Notice the extreme 
difference in the scaling of x- and y-axis in Fig. 1 relative 



to measurement accuracy: Sp » 2 x 10 m, 



10" 4 rad.) 



Due to this sensitive dependence of the rotation on the 
parameter p, a conventional experimental parameter study 
would also miss a significant part of the upper stable part 
of the branch in Fig. 1. 

In [Sieber and Krauskopf (2008)] we proposed a method 
for continuing unstable periodic orbits and bifurcations 
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Fig. 2. Pseudo block diagram for experiment with a 
tunable system parameter p and a feedback loop: 
(a) standard feedback loop PD (proportional-plus- 
derivative) with periodic reference signal 4>{i), (b) 
with projected time-delay (PTDF) block and constant 

scalar reference 4>q. The square blocks with label r 
delay the signal by t = T = Itt/lu. 

for experimental set-ups of the form shown in Fig. 2(a). 
The method assumes that a feedback loop is implemented 
around the experiment (shown as a PD (proportional- 
plus-derivative) controller in Fig. 2(a)) which achieves 
stabilization in the following sense (formulated for single- 
input-single-output feedback loops): 

(1) if the reference input <p(t) is identical to a periodic 
orbit <j)*(t;po) of the uncontrolled experiment at 

parameter p then 0*(i;po) = 0(^) is an exponentially 
stable periodic orbit of the controlled experiment. (It 
is a periodic orbit of the controlled experiment because 
the input u into the PD control (see Fig. 2(a)) vanishes 
for ((>= </>*.) 

(2) For parameter values p sa p and reference signals 
<j){t) ~ (f> t {t;po) of period T the output is asymptoti- 
cally also periodic with period T and the map 

Moo ■ K x C pcr ([0,r];]R) ^ C per ([0,T];M), (6) 
defined by 

M oo (p;^(0)(*)=^ M y(t) (7) 
is a locally well-defined and smooth map. The notation 
(j) asy (t) refers to the (periodic) output of the controlled 
experiment with inputs p and <p(t) after transients 
have died out, and C per ([0, T]; R) is the space of all 
continuous real- valued periodic functions. 

We call the feedback loop stabilizing only if both conditions 
are satisfied (possibly only locally) . That a generic feedback 
loop can be made stabilizing is a consequence of generic 
feedback stabilizability of periodic orbits [Nam and Ara- 
postathis (1992)]. The periodic orbits of the uncontrolled 
experiment can be recovered as solutions of the nonlinear 
fixed point problem 

M 0o (p;<P) = <P, (8) 
which, after discretization of (say into its N first Fourier 
coefficients), has one more variable (N + 1) than equations 
(N) and can be solved by a Newton iteration embedded 
into pseudo-arclength continuation. The Newton iteration 
defines a sequence of scalar constant inputs p and periodic 



inputs </>(•), and the residual required by the Newton 
iteration is the (periodic) asymptotic limit of the control 
signal u(-) (see Fig. 2(a)). Since (8) is an equation in the 
infinite-dimensional space R x C por the accuracy of the 
result (that is, how small ||w||oo can be made and, hence, 
how close <f> is to the unknown periodic orbit) depends on 
the choice of N. (A larger N means better approximation 
but also that a larger number of repeated experiments with 
small input variations is necessary.) 

An alternative to the Newton iteration in the space R x C per 
is shown in Fig. 2(b). The difference to the setup shown in 
Fig. 2(a) is the presence of an extra block, which we called 
"PTDF" (for projected time-delayed feedback) in Fig. 2(b), 
and which is inserted into the feedback loop (inside the 
dashed rectangle). This block implements the recursion 

4>(t) = (l-R)4(t-T)+R-(<f>(t-T)-Nrg[<f>(t-T--)]) (9) 
where < R < 1, and 

avg[#t- •)] = ^ T <Kt-8)ds (10) 

is the average of the signal <j> over the past forcing period 
T. The block feeds <j>{t) — <j){t) back into the feedback loop. 
Thus, the input u into the PD controller is 

u(t) = cj>(t) - fa) - 0o. (11) 

This is a projected version of the extended time-delayed 
feedback control (ETDFC; [Pyragas (1992); Gauthier et al. 
(1994)]): in the extreme case R = 1 it feeds back the 
difference between the output signal <j>{t) and (p(t — T) — 
avg[0(i — T — •)] + <fi , which is the output from one period 
ago but its average is shifted to the fixed input <fi . For 
R < 1 the function <j> is a weighted sum of the outputs 
from past periods (see Gauthier et al. (1994) for details). 

The feedback controlled system in Fig. 2(b) has the 
property that, whenever the output of the experiment 
is periodic (with the same period T as the delay inside 
the block "PTDF" and the forcing), the output avg[</>] is 
constant. Furthermore, if the controlled system in Fig. 2(b) 
converges to a stable periodic motion <j) c (t) with period T 
and the limit of avg[0] is identical to its scalar input 0o : 

lim a,vg[4>] = avg[0 c ] = <j> (12) 

t^oo 

then 4> c (t) is a periodic orbit of the uncontrolled experiment. 
The control signal u converges to zero if (12) is satisfied 
due to (9) and (11). 

This implies that, if the projected ETDFC system shown 
in Fig. 2(b) has a two-parameter family of stable periodic 
orbits in the (p, 0o)-planc of input parameters, we can 
define the smooth map Mi : R 2 i— » R by 

M 1 {p,4> )= lim avg[0(t--)], (13) 

t— »oo 

and, whenever the input parameters satisfy the condition 

M 1 (p;4>o) = 4>o (14) 
then the stable periodic orbit of the controlled system is 
identical to a periodic orbit of the uncontrolled experiment. 
The curve shown in Fig. 1 has been obtained as the curve 
of points in the (p; o )-plane satisfying (14). 

We note that (14) is a scalar nonlinear equation in contrast 
to the infinite-dimensional problem of the original algorithm 



proposed in [Sicber and Krauskopf (2008)]. Here Mi is 
obtained by setting the input parameters in the controlled 
system in Fig. 2(b), waiting for transients to decay and 
then measuring the asymptotic value of the scalar output 
avg [0], calling it Mi(p 7 <p ) and inserting it into (14). 

Thus, the system shown in Fig. 2(b) reduces the infinite- 
dimensional nonlinear problem, as posed by the system 
in Fig. 2(a), to a scalar equation (which gives a two- 
dimensional system after including the pseudo-arclength 
extension) at the cost of the additional block in the feedback 
loop which has to be evaluated in real-time in parallel to 
the experiment. 

5. STABILITY OF THE CONTROLLED SYSTEM 

The remaining open questions are: suppose that the 
uncontrolled experiment has a family of periodic orbits 
4>*(t;p) and assume that the feedback loop as shown in 
Fig. 2(a) is stabilizing for <p{t) = (f>*(t;p). 

(1) When is the periodic orbit (f>*(t;p) of the correspond- 
ing projected time-delayed feedback controlled system 
as shown in Fig. 2(b) also stable for </> = avg[0*]? 

(2) What is the condition of the Jacobian of the reduced 
nonlinear problem? Sometimes, reducing the dimen- 
sion of a nonlinear problem can cause a dramatic 
increase of its condition. (For example, if one reduces a 
periodic boundary- value problem to its corresponding 
fixed-point problem of the stroboscopic map.) 

Question 1 can be answered for a generalized version of the 
projected time-delay block shown in Fig. 2(b). Define the 
(Fourier) spectral projections P N : C([0,T];R) ^ R 2N+1 
and Q N : R 2N+1 ^ C([0,T];R): 

[PNV{-)]k =\\ T b k (2irs/T)y(s) ds, k = —N . . .N 
1 Jo 

N 

[Q N x]{t)= x k b k {2TTt/T) 

k=-N 

where 

b k = V2/fcos(kt) for k < 0, 

b k = V2/fsm(kt) for k > 0, and 

b Q = Vl/f. 

The set-up in Fig. 2(b) has to be generalized such that the 
input 0o is not a constant but a periodic signal given by 
the combination of the harmonic oscillators of frequencies 
up to N/T corresponding to the vector x £ R 2N+1 : 

Mt) = [QNx](t). 

Furthermore, the projected time-delay block feeds 
4>(t) = (1 - R)4>(t -T) + R[<P(t - T) - P N {<P{t — T — •))] 
back into the feedback loop and gives P N [cf)(t — •)] as its 
2N + 1-dimensional output (instead of the scalar avg[0]). 

For this general projected time-delayed feedback control 
system the following holds: the periodic orbit <j>*(t;p) of 
the uncontrolled experiment is stable in the controlled 
system for sufficiently large N, sufficiently small R, and 
input parameters <f>o(t) sufficiently close to QnPn4>*- This 
generalized time-delayed feedback may be expensive (and 
impossible to perform in real-time) if the necessary N 
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Fig. 3. Stability (a) of the projected time-delayed feedback 
system shown in Fig. 2(b) and condition of the re- 
sulting Jacobian matrix (b) for the vertically excited 
pendulum. The parameters in (3) are chosen to match 
the experimental data in [Sieber et al. (2008); see 
caption of Fig. 1]. 

is very large. It may also converge rather slowly if the 
necessary R has to be chosen close to zero. 

However, for a given experiment such as the rotations of the 
vertically excited pendulum studied in [Sieber et al. (2008)] 
the necessary N may be very small and the permissible 
R close to 1. In fact, it turns out that for the study of 
rotations near the saddle-node bifurcation N can be chosen 
equal to zero, which corresponds to the set-up in Fig. 2(b). 
Moreover, the relaxation parameter R can be chosen equal 
to 1 which simplifies the input u in the feedback loop to 

u(t) = 0(i) - <j>(t ~T)~ avg[0(i — T — •)]. (15) 

Figure 3(a) shows the stability of the periodic orbit of the 
uncontrolled system (3) when one applies the feedback loop 
shown in Fig. 2(b) with R=l. The PD box in Fig. 2 has 
input u and output 

PD[u] = -mlG-[u + 0.5u], (16) 

which, in the simulation obtaining Fig. 3, is added to 
the right-hand-side of the model (3) such that the overall 
controlled system in our simulations is 



ml 2 (f> + bcj) + bui + ml[g + uj 2 psin(ujt)]sm((j) + cut) = PD[u] 

(17) 

where PD[u] is defined by (10), (15) and (16). This is 
an idealization of the experimental set-up in [Sieber et al. 
(2008)] where control had to be superimposed with the 
up-and-down excitation. The x-axis in Fig. 3 shows the 
common factor G of the control gains in the PD control 
(16). The y-axis is the phase avg[0 + ] of the periodic rotation 
of the uncontrolled system (3). As Fig. 1 shows, near the 
saddle-node bifurcation the family of periodic rotations 
of the uncontrolled pendulum (3) cannot be parametrized 
by the system parameter p but by its phase avg[0*]. Each 
point in the plane in Fig. 3(a) shows the largest Floquet 
multiplier of the periodic orbit 0*(-;p) as a periodic orbit of 
the controlled system ( 17) . The vertical line at G = shows 
the stability of the periodic orbit without control: at G = 
and avg [</>*] = the dominant Floquet multiplier passes 
through 1. For varying G this is a transcritical bifurcation. 
Increasing gains shift this loss of stability (black curve) 
along the originally unstable part of the branch toward 
lower avg [</>*]. Figure 3 shows that the controlled system is 
stable for sufficiently large gains G. The dominant Floquet 
multiplier does not decrease uniformly for increasing G 
because we chose the ratio between proportional and 
derivative term in the PD control (16) fixed (at 0.5) and, 
thus, non-optimal. The stability chart looks similar for 
classical PD control using u — <f> — 4>„ (that is, assuming 
that we knew the periodic rotations 0* of the uncontrolled 
pendulum perfectly). 

The continuation procedure used in [Sieber et al. (2008)] to 
obtain the family of rotations in Fig. 1 starts from a stable 
rotation at p ss 2 cm where one one can simply measure 
avg [(/>*] and assign the initial input 0o to this value. The 
initial unit tangent is (ptan, 0o,tan) = ( — 1,0). Then in each 
continuation step one performs a Newton iteration, running 
a sequence of controlled experiments as shown in Fig. 2(b) 
for a sequence of inputs (p, 0o) as required by the Newton 
iteration and measuring the residual r = (r%, r 2 ) given by: 

r% = Ptan ' (P ~ Paid) + 0O,tan ' (00 ~ 0O,old) - h 

r 2 = M 1 (p; O ) - O 
where (p ld, 00 old) is the point previously found in the 
continuation, and (ptan, 0o,tan) is the (approximate) unit 
tangent to the solution curve in (p id, 0o,oid)- The Newton 
iteration is successful if the norm of the residual is smaller 
than a given tolerance (5 x 10~ 3 in [Sieber et al. (2008)]). 
The evaluation of r 2 requires running the controlled 
experiment until the transients have decayed. (How long 
this takes can be estimated from Fig. 3(a).) 

This leads to the second question: how robustly does the 
Newton iteration converge? The convergence of the Newton 
iteration depends on the condition cond(J) = || J|| ■ || J _1 | 
of the Jacobian J = [dr /dp, dr / d<fio) , which is shown in 
Fig. 3(b). We note that || J~ 1 \\ is always 1 for the system 
parameters in [Sieber et al. (2008)] (or very close to unity 
if the tangent is only approximate) because the rows of J 
are orthogonal to each other by definition of the tangent 
(ptan, 0o, tan), and dr 2 /dp is large (due to the sensitive 
dependence of the rotation on the excitation amplitude 
and the moderate decay rate of the controlled system as 
shown in Fig. 3(a)). Figure 3 provides evidence that the 



nonlinear system is uniformly well-conditioned for a wide 
range of gains. The only effect increasing the condition of 
the Jacobian is the loss of control when the gain becomes 
too small near the transcritical bifurcation (black curve in 
Fig. 3(a) and (b)). 

6. CONCLUSION 

We have analysed the stability and robustness of the 
experimental continuation of the periodic rotations of a 
pendulum through a saddle-node bifurcation performed in 
[Sieber et al. (2008)]. This analysis is important because 
the experiment relies on the asymptotic convergence of 
a feedback controlled experiment and Newton iterations, 
and in a real experiment one can never achieve that the 
input u into the PD controller, defined by equation (11), 
vanishes perfectly (due to disturbances and incomplete 
decay of transients, and a non-zero tolerance of the Newton 
iteration) . We perform our analysis for a model equation. 
Because one is interested in the order of magnitude of 
condition numbers and decay rates. This approach is 
justified as long as the model is qualitatively correct (which 
is the case for a menchanical pendulum) . We showed that 
both parts of the continuation process — the experiment 
with its projected time-delayed feedback control and the 
Newton iteration in K 2 for the inputs into the controlled 
experiment — converge uniformly in the vicinity of the 
saddle-node bifurcation and on the unstable part of the 
family of rotations. 
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